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Background: 50 % of the heavy element abundances are produced via slow neutron capture 
reactions in different stellar scenarios. The underlying nucleosynthesis models need the input of 
neutron capture cross sections. 

Purpose: One of the fundamental signatures for active nucleosynthesis in our galaxy is the obser¬ 
vation of long-lived radioactive isotopes, such as ®°Fe with a half-life of 2.60 x 10® yr. To reproduce 
this 7-activity in the universe, the nucleosynthesis of ®®Fe has to be understood reliably. 

Methods: A ®®Fe sample produced at the Paul-Scherrer-Institut was activated with thermal and 
epithermal neutrons at the research reactor at the Johannes Gutenberg-Universitat Mainz. 

Results: The thermal neutron capture cross section has been measured for the first time to 
(Jth = 0.226 (1® Q49) b. An upper limit of ctri < 0.50 b could be determined for the resonance 
integral. 

Conclusions: An extrapolation towards the astrophysicaly interesting energy regime between 
kT = 10 keV and 100 keV illustrates that the s-wave part of the direct capture component can 
be neglected. 

PACS numbers: 25.40.Lw, 26.20.Kn, 27.50.+e, 28.20.Ka 


I. INTRODUCTION 

The decays of the unstable isotopes ®°Fe 
(ti /2 = 2.60Myr [T]) and ^®A1 {ti /2 = 0.717Myr [5]) 
in the Milky Way, which have been observed with 
satellite-based 7 -ray telescopes 011], are considered as 
a clear signature of ongoing stellar nucleosynthesis 0 . 

The production of ®°Fe in the slow neutron capture 
process (s-process) 0 is hampered by the rather short¬ 
lived precursor ®®Fe {ti /2 = 44.495 d 0), which acts as a 
branch point of the s-process path as illustrated in Fig.[l] 
Accordingly, high neutron densities are required to avoid 
that the reaction flow is bypasses ®°Fe via the decay of 
®®Fe. Once ®°Fe is reached, it can also be destroyed by 
neutron capture or - on longer time scales - by /3“-decay. 
High neutron densities are generally accompanied by very 
high temperatures, but the synthesis of ®°Fe requires an 
upper limit of about 2 x 10® K (Tg = 2), because photo¬ 
disintegration reactions such as ®®Fe( 7 , n) and ®®Fe( 7 , n) 
start to dominate otherwise. 

There are two different astrophysical scenarios where 

60 pe 

can be produced 0: during the He-shell burning 
phase in low-mass thermally pulsing asymptotic giant 
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FIG. 1: (Color online) The s-process reaction path between 
Fe and Ni. The isotope ®®Fe is produced via a sequence of 
(n,7) reactions starting at the stable iron isotopes. Because 
of the short half-life of ®®Fe (ti/2 ~ 45 d), the production of 
®®Fe depends critically on the stellar neutron density. 

branch (AGB) stars and during the convective C-shell 
burning in massive pre-supernova stars. In AGB stars, 
neutron densities of 10 ^® cm“^ and temperatures around 
2.5 X 10® K (Tg = 2.5) are reached, whereas in mas¬ 
sive stars neutron densities of 10 ^® cm“® at temperatures 
of up to Tg = 10 during G-shell burning are reached 
[10] . According to detailed stellar model calculations by 
Limongi and GhiefR 0 , about 65 % of the total yield of 
®®Fe are in fact synthesized in the pre-supernova stage of 
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FIG. 2: (Color online) The detection efficiency of the HPGe detector used for the activity measurements in Mainz. The 
efficiency used for determination the number of ®°Fe particles is given in Table H] The solid line shows the least squares fit 
to interpolate between the data points of the calibrated solution (red). For and ®°Co the data points were corrected for 
cascade effects using the GEANT-3.21 package [Si]. The grey band represents the uncertainty of the fit. 


massive stars and 18 % are contributed by the He burn¬ 
ing shell of less massive stars. A third major component 
is eventually produced by explosive shell burning dur¬ 
ing the supernova itself. These contributions to the total 
®°Fe yield are strongly affected by the respective masses 
and metallicities of the stars involved and may vary cor¬ 
respondingly. 

A crucial input for the production of ®°Fe in AGB stars 
and massive pre-supernova stars are the neutron cap¬ 
ture cross sections at the respective stellar temperatures. 
So far, an activation measurement of the ®°Fe(n, 7 )®^Fe 
cross section at neutron energies corresponding to a ther¬ 
mal energy of kT = 25keV (typical for AGB stars) 
was performed at Forschungszentrum Karlsruhe, Ger¬ 
many. The Maxwellian averaged cross section (MAGS) 
at kT = SOkeV was determined to (5.15 ± 1.4) mb [TT] . 
The direct capture (DG) component of the cross section 
at this temperature constitutes an important information 
for the extrapolation towards the astrophysically inter¬ 
esting temperatures in massive stars around kT = 90 keV. 
In this respect, the thermal cross section provides a con¬ 
straint for the s-wave component of the DG cross section. 
Therefore, the previously unknown thermal cross section 
of ^'^Fe was measured using the irradiation facility at the 
TRIGA (Training, Research, Isotopes, General Atomic) 
type research reactor at Johannes Gutenberg-Universitat 
Mainz, Germany [nisi. 

II. EXPERIMENT 

The ®°Fe sample was produced at the Paul-Scherrer- 
Insitut (PSI) in Villigen, Switzerland [14]. In order to 
compensate for the limited amount of ®°Fe the only pos¬ 
sible method for the determination of the neutron cap¬ 


ture cross section was an integral activation measure¬ 
ment at high neutron fluxes. Gompared to the more 
generally applicable time-of-flight technique, the activa¬ 
tion method has the advantage of excellent sensitivity 
m, which allows neutron capture measurements even 
on very small samples [m ds). This technique has 
the additional advantage that it does not require iso- 
topically enriched samples, because the capture reac¬ 
tions can be identihed via the 7 -decay characteristics 
of the product nucleus ®^Fe. 

A. Measurements and calibration 

The induced activities were measured using a HPGe 
detector (GANBERRA-GX7020) with a relative effi¬ 
ciency of 72.3%. The output signals from the pream¬ 
plifier were converted with a flash-ADG (CAEN module 
V1724). The dead time corrections were determined us¬ 
ing a ^^^Cs sample, which was placed at a fixed distant 
position during all activity measurements. The corre¬ 
sponding corrections were negligibly small. Because of 
the contamination of ^®Fe in the ®°Fe sample, the activity 
of ^®Fe was suppressed by a lead foil 1 mm in thickness. 

The efficiency was determined by a calibrated solu¬ 
tion containing the standard single- or double-line 7 -ray 
emitters ® 0 Co, ^^Sr, 8 ®Y, and ^osHg. 

The uncertainty of the 7 -emission rate was given with 3 % 
(2 (t). This multi-nuclide solution was absorbed in a pure 
graphite disc 6 mm in diameter and 1 mm in thickness to 
match the properties of the ®°Fe sample used in the mea¬ 
surement (see below). For all 7 -activity measurements, 
the samples were placed 7.4 mm in front of the Ge crystal. 
Because of the small distance between sample and detec¬ 
tor, cascade corrections were necessary for the decays 
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TABLE I: Decay characteristics and efficiency of 7 -emission of the investigated nuclei. 


Isotope 

tl/2 

E., /keV 

c /% 

Reference 


®°Co 

(1925.28 ±0.14) d 

1173.228 ± 0.003 

99.85 ±0.03 

[13 

0.024 ±0.0002“ 



1332.492 ± 0.004 

99.9826 ±0.0006 


0.022 ± 0 . 0002 “ 

^^Zr 

(16.749 ±0.008) h 

743.36 ± 0.03 

93.09 ±0.16 

m 

0.069 ± 0.002® 

®®Zr 

(64.032 ± 0.006) d 

724.195 ±0.004 

44.27 ±0.22 

m 

0.070 ±0.001® 



756.728 ±0.012 

54.38 ±0.22 


0.068 ± 0 . 002 ® 

®iFe 

(5.98 ± 0.06) min 

297.90 ±0.07 

22.24 ± 2.88 

m 

0.087 ±0.002® 



1027.42 ±0.11 

42.73 ± 4.92 


0.056 ±0.002® 



1205.07 ±0.12 

43.60 ± 4.50 


0.051 ±0.001® 


“ measured with a HPGe detector at Goethe-Universitat Frankfurt (used for 60Fe determination). 

^ measured with a different HPGe detector at the research reactor at Johannes Gutenberg-Universitat Mainz (used for ®^Fe 

and Zirconium determination) 


of ®°Co and ®®Y. Those corrections were based on the 
simulations performed using the GEANT-3.21 package 
[311]. The corresponding correction for the emission line 
of ®°Co at the energy of 1173 keV was 30%, at 1332 keV 
31 % and for at the energy of 898 keV and 1836 keV 
27% and 29%, respectively. As shown in Fig.j^the mea¬ 
sured efficiencies could be reproduced within the exper¬ 
imental uncertainties of ±3.5 % over the energy range 
from 150 keV to 1900 keV by the expression 

e-y = A exp[ — B ln{i?.y — C + D exp(F x i^T-)}]. (1) 


B. Sample preparation 

The ®°Fe was extracted from slices of a cylindrical cop¬ 
per beam dump, which was previously irradiated with 
590 MeV protons at PSl [Uj. In addition to ®°Fe activity, 
the initial copper sample of 3 g also contained 150 MBq of 


0 + 


0.0 


'‘"Fe \ 

(2.60 Myr) 



FIG. 3: (Color online) The activity of the ®®Fe sample was 
determined by the 7 -ray cascade with energies of 1173 keV and 
1332 keV emitted after the /3-decay of the daughter nucleus 
®®Co. Data from [13 mi- 


®°Co, 100 MBq of ®®Fe, and 2 MBq of "'^'^Ti. Details of the 
chemical separation of the ®°Fe fraction are described in 
M- The final purification was performed shortly be¬ 
fore the experiment using liquid-liquid extraction into 
methyl-isobutyl ketone from 7 M HCl solution and fol¬ 
lowing back-extraction with diluted HCl. This solution 
was evaporated on a graphite disc with 6 mm diameter 
and 1 mm thickness. The number of ®°Fe atoms in the 
sample was determined via the increasing ^'^Co activity 
(ti /2 = 5.272y [13) according to Fig. The activity 
of ®°Co nuclei increases as 

A(®°Co) = (1 - exp[-A(®°Co)t]) A(®°Fe), (2) 

where A is the decay constant. The related 7 - 
activity at time t can be derived using the integrated 
line strength C.y 



Ey [keV] 


FIG. 4: (Color online) 7 -ray spectrum for the determination 
of the ®®Co activity taken 34 months after purification of the 
sample. The measurement time was tm = 23 h. 
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TABLE II: Number of atoms N (in units of 10^®) and cross 
sections (in barn) of the neutron fluence monitors ®'^Zr and 
®«Zr . 



®^Zr 

®®Zr 

N 

2.7331 ± 0.044 

0.440 ±0.014 

Nci 

2.7101 ± 0.044 

0.437 ±0.014 

O'th 

0.0494 ±0.0017 

0.0229 ±0.0010 

ctri 

0.280 ±0.010 

5.28 ±0.11 


Cross sections were obtained from Ref. |22[ . 


TABLE III: The amount of Zr nuclei produced in the acti¬ 
vations and the resulting neutron fiuences without and with 
cadmium shielding^. 




without 

Cd 


with Cd 


iV(95 

Zr) 

1.517 ±0.005 

±0.018 

0.510 

± 0.006 ± 

0.008 

iV(®'^ 

Zr) 

1.171 ±0.001 

±0.018 

1.067 

± 0.001 ± 

0.016 

<E>th ' 


8.60 ±0.03 

±0.38 

1.21 

±0.01±i 

0.16 

4>epi 

C 

0.467 ±0.002 

±0.014 

0.458 

± 0.005 ± 

0.014 


“ Uncertainties are statistical and systematic, respectively. 
^ In units of 10®. 

In units of 10^^cm“®. 


AtCCo) 


0.9976 

_A_ 

exp{-Al} - exp{-A(I -b t„i)} ’ 


(3) 


where A is the decay constant of ®°Co and with the 
factor 0.9976±0.0003 [T7] for the fraction of ®°Co”*, that 
decays to the ground state of ®°Co, the measurement time 
and the detection efficiencies for the 1173 keV and 
1332 keV y-transitions, respectively. For the analysis of 
the emission line at the energy of 1332 keV, the decay of 
®°Co™ has to be corrected. The decay intensities and 
efficiencies are listed in Table |Tj which summarizes all 
decay characteristics adopted in the data analysis. With 
A = XN, the number of ®°Fe atoms becomes 


7V(60Fe) 


At(^°Co) 1 

1 — exp{—A(®°Co)t} A(®°Fe) 


(4) 


The activity measurement of ®°Co was carried out at 
the GoetheUniversitat Frankfurt 38 months after the pu¬ 
rification using an HPGe detector of 98% relative effi¬ 
ciency (see Fig. |^. Background due to the activity of 
the ®®Fe contamination in the sample was suppressed by 
a lead foil 1 mm in thickness. The number of ®°Fe atoms 
in the sample 


Af(®°Fe) = (7.77 ± 0.11,,,. ± 0.42,,„) x 10^^ (5) 


has been determined as a weighted average comprising 
both ®°Go lines. The systematic uncertainty is deter¬ 
mined by the y-ray detection efficiency, the decay inten¬ 
sities, and the half-life (see Table [I]). As the half-life of 
®°Fe a value ti /2 = (2.60 ± 0.05) Myr [T] was used. 


C. Reactor activations 


The ®°Fe sample was activated for = 10 min with 
and without cadmium foils surrounding the sample in 
both cases. This so-called cadmium-difference-method 
allows the distinction between the thermal neutron cap¬ 
ture cross section and the resonance integral, which 
takes into account the epithermal component of the 
reactor neutron spectrum. The reactor spectrum can 
be described as the sum of a thermal component, 
i.e., a Maxwell-Boltzmann distribution corresponding to 
kT = 25.3 meV, and an epithermal component follow¬ 
ing an l/i?-dependence. Due to the very large thermal 
capture cross section of cadmium, a proper cadmium 
shielding of the sample results in a significantly differ¬ 
ent response to thermal and epithermal neutrons. In the 
ideal case, all thermal neutrons would be absorbed in the 
cadmium, while the epithermal spectrum remains undis¬ 
turbed. 

The number of product nuclei after the activa¬ 
tion N{^+^X) can be expressed in terms of the 
thermal cross section ath, the resonance integral 



[keV] 


In view of the short half-life of the produced ®^Fe nu¬ 
clei (ti /2 = 5.98min |20)). the activations at the TRIGA 
research reactor were performed using a pneumatic trans¬ 
port system between the irradiation position and the 
counting room [H [SI- 


FIG. 5: (Color online) The y-ray spectra of the Zr moni¬ 
tor foils used for the neutron fluence analysis normalized for 
different measuring times. Due to the smaller capture cross 
section of ®®Zr in the thermal energy regime, the cadmium 
shielding affects the ®^Zr signature to a much lower degree 
than that of ®®Zr. 

















5 


o’Ri (ctri = cr(-B)/-Bdi? with the cutoff energy 

Ecutoff ~ 90 meV), and the epithermal ($epi) and ther¬ 
mal neutron ffuences (d>th) in units of cm“^, 

N{^+^X) = N{^X) ($thath + ^epifTRi), (6) 

where N{^X) is the number of target nuclei in the irra¬ 
diated sample. Natural zirconium provides a well suited 
monitor for the epithermal and the thermal flux. The ac¬ 
tivation of the isotopes ®"^Zr and ®®Zr exhibit significantly 
different ratios crRi/dth (Table [ll|. The uncertainties in 
the number of Zr atoms are due to the sample weight 
(0.2%) and to the isotopic abundances (1.6 and 3.2% 
for ®'^Zr and ®®Zr, respectively) [53]. Two sets of Zr foils 
were used in the activations with and without cadmium 
shielding. 

III. RESULTS 


TABLE IV: The number of ®^Fe nuclei (in units of 10®) pro¬ 
duced in the activations. 


7 -ray energy 

fV(®^Fe)‘ 

1 

/keV 

without Cd 

with Cd 

1027 

1.54 ±0.19 ±0.18 

< 0.179 

1205 

1.48 ±0.20 ±0.16 

< 0.206 

Weighted average 

1.51 ±0.14 ±0.24 

< 0.179® 

“ Uncertainties are 

statistical and systematic. 

respectively. 


® Adopted upper limit for further discussion. 

of 0.5 %. The systematic uncertainty is again determined 
the 7 -efficiencies, the decay intensities, and the half-life 
(Table |T|. The resulting neutron ffuences for the two ac¬ 
tivations are listed in Table uni 


A. Determination of the neutron fluence 


B. Thermal ( 0 , 7 ) cross section of ®°Fe 


The zirconium foils used as neutron monitors are 
0.127 mm in thickness and 6 mm in diameter. The foils 
are thin enough that neutron self-absorption losses dur¬ 
ing 7 -spectroscopy could be neglected. The ffuences in 
units of 1 /cm^ are 


^epi — 


Ar97 _ 


th'i'th 




.96 

’^RI 


(7) 


and 


where the indices are referring to the various Zr iso¬ 
topes. Figure shows the 7 -ray spectra of the monitor 
foils normalized to equal neutron fluence. Because of the 
small neutron capture cross section of ®®Zr in the thermal 
energy regime, the ®^Zr signal is only marginally affected 
by the cadmium shielding, whereas ®®Zr exhibits a clear 
effect due to the larger thermal cross section of ®^Zr. The 
number of produced Zr nuclei is 


N{^X) 






where 


_ 1 - exp(-A»ta) 

/a , 

fw —exp ( Ajtw)! 

fra —1 exp ( Xityn) 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


are the corrections for the decay during the activation 
/a, during the waiting time between activation and mea¬ 
surement /w, and during the measurement /m, respec¬ 
tively. The correction for the deadtime was of the order 


The 7 -spectrum measured after the activation of the 
®°Fe sample without cadmium shielding (Fig. clearly 
exhibits the 7 -transitions of ®^Fe at 297.9keV, 1027keV, 
and 1205 keV. However, only the last two ones were used 
in the analysis because of the poor signal-to-background 
ratio of the 298 keV line. The systematic uncertainty is 
calculated by the error of the efficiency, the Ij, the half- 
lifes, and the neutron ffuences (see TablejTj Table | 111 [ and 
Table [i^. 

In the corresponding spectrum measured after the acti¬ 
vation with cadmium shielding, the ®^Fe lines are com¬ 
pletely missing as illustrated in Fig. for the 1027 keV 
line as an example. In this case, only an upper limit can 
be determined for the resonance integral. Th e numbers 


of produced °^Fe nuclei are listed in Table IV 


The number ratio of ®^Fe and ®°Fe after the activation 
without cadmium is 


7V(6iFe)/A(60Fe) = + ^epi^Ri- 


The thermal cross section 


ath(®°Fe) 


iV(6lFe) 1 $epi 

iV(60Fe) $th $th 


( 12 ) 


(13) 


is determined by the number of sample atoms iV(®°Fe) 
(Sec. 11B), the neutron ffuences $th and $epi from the 


Zr monitor measurements (Table 111), and the number 
of ®^Fe nuclei produced during the activations N{^^Fe) 


(Table [rv| . The resonance integral 


N'- 


He) 


ctri = 


Af(f*0Fe) 


SS' NC^Fe) 

$th N{B0Fe) 


(fjCd _ ‘t’Sf *^6?: 

epi 


(14) 


is obtained accordingly. Since the epithermal ffuences 
were almost equal in both activations, and because the 
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FIG. 6 : (Color online) The 7 -ray spectrum of the activation without cadmium shielding measured for one, two, and three 
half-lives of ®^Fe, respectively. The upper panel provides an overview, and the lower panels show a zoom into the regions of 
the ®^Fe lines at 1027keV and 1205 keV. The background stems from the activation of the contaminants of the ®°Fe sample. 


number of nuclei produced with the cadmium ab¬ 
sorber is much smaller than without absorber, an upper 
limit for the resonance integral can be defined as 


^Cd( 61 pg) 2 

jV( 60 Fe) 


(15) 


Assuming a Icr confidence level as a constraint for the 


resonance integral derived from the 1027 keV line, one 
finds 


for calculating the thermal cross section using Eq. (131. 


A variation of the resonance integral within these lim¬ 
its affects the thermal cross section by about 10%. We 
assume the resonance integral here explicitly as 



I 7x10- 

^ 5x10- 


970 


1027 

[keV] 


1100 


am = 0.00l°;^° b, 


(17) 


consistent with Eq. and treat this range as a sys¬ 

tematic uncertainty. Should the resonance integral be 
improved in the future, the thermal cross section can be 
re-evaluated accordingly. With Eqs. (131 and 0 the 
thermal cross section of ®°Fe becomes 

ath(®°Fe) = (0.226 ± 0.021,tat(to;o4l)syst) b. (18) 


IV. SUMMARY AND DISCUSSION 


FIG. 7: (Golor online) A detailed view into the region around 
1027 keV of the 7 -ray spectrum after activation with the cad¬ 
mium shielding, illustrating the absence of ®^Fe lines. There¬ 
fore, only an upper limit can be deduced from the data. 


Within this work, we characterized the ®°Fe sample to 
contain A(®°Fe) = (7.77±0.11„tat±0.42^y„t) x 10^^ atoms. 
Using the cadmium-difference-method two activations of 
that sample have been performed at the TRIGA research 
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FIG. 8: (Color online) The Maxwellian averaged cross section for ®°Fe(n, 7 ). The present measurement of the cross section 
at kT = 25.3 meV (red triangle) and an 0.00115/\/S extrapolation to the astrophysical energy regime are indicated by the 
solid red triangle and the dashed red line, respectively. This extrapolation can be used to estimate the DC of the MACS at 
kT = 25keV (black dot) [11]. The astrophysical energy regime from kT = lOkeV - 100 keV (grey box) is clearly dominated 
by the resonant capture contribution. Below about 1 keV, the MACS based on the most recent version of the TENDL library 
(TENDL-2014 |24| . blue line) are a factor of 3 above the current measurement. This indicates that the DC componenent is 
clearly overestimated in the library. 


reactor at Johannes Gutenberg-Universitat Mainz, Ger¬ 
many. The neutron capture cross section of ®°Fe at ther¬ 
mal energies and an experimental upper limit for the res¬ 
onance integral could be determined for the first time: 

ath(®°Fe) = 0.226 (^^ 049 ) b 

and 

ctri < 0.50 b. 

Figurej^shows a comparison of our data with evaluated 
cross sections (TENDL-2014 [24]) and the so far only ex¬ 
perimental value of 5.7±1.4 mb at fcT = 25 keV [mi25j . 
Under the assumption that the MACS in the meV-regime 
is dominated by an s-wave direct capture component, an 
extrapolation towards higher energies via ljV~E is possi¬ 
ble. Together with the measurement of the total capture 
cross section at kT= 25keV, it is then possible to dis¬ 
entangle the direct and the resonant contribution in the 


astrophysically interesting energy regime. It turns out 
that the direct capture component is almost negligible, 
ranging from less than 10 % to less than 1 % between 
lOkeV and 100 keV. The comparison of the experimen¬ 
tal data with the latest release of TENDL indicates that 
the resonant component is well described, but the direct 
capture component is overestimated. 
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